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ABSTRACT 

The  scanning  tunneling  microscope  (STM)  was  used  to  quantitatively  investigate 
sputter  etching  of  graphite  surfaces  bombarded  by  5  keV  argon  ions.  The 
resulting  surface  morphology  depended  strongly  on  ion  flux,  dose  and  sample 
temperature.  The  height-height  correlation  functions  of  the  roughened  surfaces 
were  calculated  directly  from  the  STM  topographs  and  were  compared  to  linear 
response  theory  and  scaling  analyses  for  the  propagation  of  growing  fronts.  We 
find  that  the  surfaces  develop  structures  characterized  by  a  correlation  length 
that  diverges  with  increasing  ion  dose. 
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One  of  the  fundamental  challenges  in  materials  science  is  to  understand 
the  effects  of  particle  radiation  on  solid  surfaces.  For  example,  the  properties 
of  sputter-deposited  thin  films  depend  strongly  on  the  surface  topology,  which  in 
turn  is  determined  by  the  sputter-growth  process.^  Ion-beam  sputter  etching,  a 
widely  used  technique  in  surface  science  and  microelectronics,  also  produces 
surfaces  with  characteristic  morphologies  that  depend  on  the  sputtering 
conditions. 

On  the  macroscopic  level  (i.e. ,  length  scales  larger  than  one  micron),  the 

effects  of  particle  radiation  on  surfaces  are  well  understood.  Ion  bombardment 

erodes  edges,  corners  and  peaks  preferentially  because  prominent  features  are 

exposed  to  a  larger  flux.  This  smoothing  of  surface  roughness  is  further 

enhanced  by  diffusion  and  evaporation/recondensation.  In  the  continuum 

description  of  this  process,  one  either  assumes  that  the  local  erosion  rates  are 

proportional  to  the  local  incoming  flux,  or  one  uses  the  so-called  "Huyghens 

Construction."®  Numerical  studies  of  macroscopic  evolution  have  obtained 

fi  7 

results  that  compare  favorably  with  experiments.  ’ 

At  the  sub-micron  level,  our  understanding  is  much  less  complete.® 
Electron-microscope  studies  indicate  that  particle  radiation  roughens  a  surface, 
and  prominent  cone-like  features  have  been  reported.®’*®  Both  analytical  studies 
and  Monte  Carlo  simulations  1 L12  of  the  effect  of  shot  noise  on  continuum 
theories  indicate  that  stochastic  processes  roughen  surfaces.  The  roughening  is 
characterized  by  a  correlation  length  £(t)  that  diverges  with  etching  time  t.  Until 
now,  however,  there  has  been  little  experimental  confirmation  of  this  description. 

In  this  Letter,  we  report  on  our  use  of  a  Scanning  Tunneling  Microscope 
(STM)*3  to  study  the  surface  topography  of  graphite  after  sputter  etching.  The 
STM  offers  unique  opportunities  for  the  study  of  radiation  erosion  because,  unlike 
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the  electron  microscope,  it  can  quantitatively  measure  the  height  profile  h(r)  of 
the  surface  over  a  wide  range  of  length  scales.  This  allows  us  to  calculate  the 
height-height  correlation  function^  ( |  h(q)  1 2 )  t  defined  as: 


<lh(q)|2> 


eiqr  <(h(o)-h(r))2)t/Area, 


(1) 


with  (•••  )t  indicating  a  sample  average  after  t  seconds  of  exposure.  In  other 
roughening  problems,  e.g.,  thermal  roughening,  the  corresponding  correlation 
function  is  a  good  measure  of  the  overall  surface  geometry.  15  Because  of  loss  of 
phase  information  in  the  height-height  correlation  function,  however,  it  is  less 
sensitive  to  uncorrelated  -  but  prominent  --  surface  features,  which  are  much 
more  easily  visible  in  the  real  space  STM  topographs. 

The  surface  chosen  for  this  investigation  was  the  cleaved  (0001)  face  of 
highly  oriented  pyrolytic  graphite  (HOPG).  This  surface  is  inert  in  air  and  is 
easily  imaged  with  the  STM. Graphite  also  has  a  rigid  lattice,  with  a  melting 
temperature  of  ~3800°C.  This  indicates  that  surface  diffusion  should  be  minimal 
at  room  temperature  and  thus  bombardment  induced  topography  is  "frozen  in" 
and  can  be  observed  with  the  STM  long  after  sputtering  has  occurred. 

Freshly  cleaved  graphite  samples  were  examined  with  the  STM  before 
sputter  etching.  The  microscope  was  operated  at  atmospheric  pressure  in  the 
constant  current  mode,  with  a  tunnel  current  of  0.5  nA  and  a  sample  to  tip  bias  of 
•100  mV.  No  filtering  or  data  enhancement  was  necessary  because  of  the 
especially  high  stability  of  our  STM  with  respect  to  vibrational,  electrical  and 
thermal  disturbances.  Low  magnification  topographs  (2400  A  x  2400  A  image  size) 
showed  large,  atomically  flat  areas  over  many  thousands  of  square  Angstroms, 
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while  at  higher  magnifications  (25  A  x  25  A)  the  familiar  atomic  scale  features  of 
dean  graphite  were  easily  observed  A3 

After  stable  images  of  clean  graphite  were  obtained,  the  samples  were 
transferred  to  the  sample  treatment  chamber  of  a  KRATOS  XSAM-800  surface 
analytical  system.  The  graphite  surfaces  were  sputter  etched  with  a  beam  of  5 
keV  Ar+  ions,  rastered  over  a  9  mm;  area  on  the  sample  and  inrident  at  an  angle 
of  60°  to  the  surface  normal.  The  beam  flux  incident  on  the  sample  was 
determined  by  using  an  electrometer  to  measure  the  ion  beam  current.  A  small 
positive  bias  (45  volts)  was  applied  to  the  sample  to  suppress  secondary  electron 
emission.  The  experimental  parameters  that  were  varied  in  this  study  were  the 
flux  J,  the  ion  dose  Q  =  Jt,  and  the  substrate  temperature  T. 

This  Letter  highlights  the  results  obtained  from  over  1,000  STM  topographs 
collected  from  more  than  15  bombarded  samples.  The  two  fluxes  reported  here 
are  6.9  x  10*3  ions/cm2  sec  and  3.5  x  10 14  ions/cm^  sec.  By  varying  the  time  of 
exposure  to  the  ion  beam,  the  total  doses  obtained  were  lO1^  iol7  and  101® 
ions/cm^.  In  addition  to  the  ambient  temperature  experiments,  etching  was  also 
performed  at  surface  temperatures  of  approximately  600K  and  900K  for  a  flux  of 
3.5  x  101*  ions/cm2  sec  and  a  dose  of  10*7  ions/cm2. 

The  graphite  samples  were  re-examined  with  the  STM  after  etching  using 
identical  operating  parameters  and,  if  possible,  the  same  tunneling  tip  used  prior 
to  sputtering.  The  results  shown  in  this  study  were  reproducible  from  sample  to 
sample  and  even  with  different  tunneling  tips,  demonstrating  that  tip  artifacts 
were  minimal.  Topographs  with  obvious  imaging  artifacts  were  observed,  but 
were  discarded.  Figs.  1-3  show  2400A  x  2400A  images  of  sputtered  graphite  for 
the  three  different  ion  doses.  The  sample  was  etched  at  the  lower  flux.  The 
corresponding  correlation  functions  are  shown  in  Fig.  4a.  The  sample  average 
( I  h(q)  I  )  of  the  two-dimensional  Fourier  transform  power  spectrum  was  obtained 
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by  performing  a  rotational  average  of  I  h(q)  1 2,  sampling  at  increments  of  1°  for 
each  value  of  q.  The  variance  of  <  I  h(q)  I  was  also  calculated  to  provide  an 
estimate  of  the  uncertainty  in  the  correlation  function. 

The  surface  topology  depended  strongly  on  the  ion  dose.  If  we  define  the 
corrugation  as  the  slope  of  a  line  connecting  two  surface  points  separated  laterally 
by  a  distance  L,  Fig.  1  shows  that  this  quantity  becomes  small  for  L  larger  than 
the  correlation  length,  £  «  100  •  200A.  Thus,  for  such  length  scales,  the  surface 
can  be  considered  to  be  flat.  As  the  dose  Q  increases  (Figs.  2  &  3),  %  increases  as 
well,  and  at  the  highest  doses  E,  exceeds  the  image  size.  For  this  real  space 
behavior,  one  would  expect  that  the  correlation  function  should  be  q-independent 
for  q  <  l/£,  while  it  should  decrease  with  q  for  q  >  1 /£,  as  observed  in  Fig.  4a. 

The  creation  of  correlated  structures  by  particle  radiation  is  counter¬ 
intuitive.  However,  theoretical  studies  of  the  non-equilibrium  growth  of 
interfaces1 1*12  indicate  that  correlation  resulting  from  random  events  is  possible. 
A  linear  response  theory  was  proposed  for  radiation  erosion  in  ref.  16  and  we  now 
test  its  predictions  by  comparing  the  asymptotic  limits  to  the  experimental  data. 

In  linear  response  theory: 

dth(q,t)  =  -  co(q)  h(q,t)  +  Tj(q,t)  ,  (2) 

where  co(q)  is  the  healing  rate  of  a  surface  modulation  of  wavevector  q  and  T|(q,t) 
is  the  Gaussian  white  noise  for  the  incident  ions  with  a  variance  proportional  to 
the  flux  J.  For  the  case  of  isotropic  radiation  erosion  including  annealing  by 
surface  diffusion,  the  healing  rate  oXq) «  J I  q  I  +  D I  q  1 4  ,  with  D  proportional  to 
the  surface  diffusion  constant.  The  resulting  correlation  function  is  : 

<  I  h(q)  1 2  )t  «  JL  [i .  exp  (-  2co(q)t)] .  (3) 
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According  to  Eq.(3),  for  small  q,  (  I  h(q)  I  2)t  is  proportional  to  t  and  independent  of 
q,  while  for  large  q  it  decreases  with  q,  indicating  two  distinct  regions  in  plots  of 
(  I  h(q)  I  2>t  vs.  q.  Fig.  4a  shows  that  this  type  of  behavior  is  observed 
experimentally.  The  transition  should  occur  at  the  crossover  wavevector  q0  = 
defined  by  co(£‘l)  t  =  1,  which  predicts  that  the  correlation  length  at  crossover,  £  ~ 

Q  =  Jt  (assuming  E,  <  (D/J)^  ).  From  Fig.  4a,  we  see  that  %  increases  with  dose  Q, 
albeit  more  slowly  than  linearly.  In  Fig.  4b,  we  show  ( I  h(q)  I  2>t  for  the  same 
three  doses  as  for  Fig.  4a,  but  with  J  increased  by  a  factor  of  5.  Indeed,  within  the 
uncertainties  of  the  experimental  data,  %  does  not  appear  to  have  changed 
significantly  for  surfaces  sputtered  with  the  higher  flux,  even  though  these 
surfaces  were  rougher,  ie.  had  larger  values  of  ( I  h(q)  I  2)t. 

We  also  investigated  the  behavior  of  Eq.(3)  for  elevated  temperatures.  For 
large  q,  <o(q) «=  Dq4,  so  that  I  h(q) 1 2>t  «  J/(Dq4).  This  q*4  tail  should  be  more 
pronounced  at  elevated  temperatures,  since  D  increases  rapidly  with  T.  In  Fig. 
4c,  we  show  the  T-dependence  of  <  I  h(q)  1 2>  at  the  values  of  J  and  Q  of  Fig.  2. 

n 

Above  600°K,  ( I  h(q)  I  )  drops  more  sharply  with  q  and,  for  large  q,  has  a  tail  with 
an  approximately  q'4  dependence.  For  lower  T,  we  found  no  q"4  tail  down  to  0.1 
A'l  in  any  of  our  topographs. 

By  expanding  the  exponential  term  in  Eq.(3)  for  the  limit  q  0,  we  see  that 
the  interface  width,  W  *  qS%)(  I  h(q)  1 2)tl/2  t  should  be  proportional  to  (Jt)^2  =  Ql/2( 
a  prediction  which  is  independent  of  our  choice  of  the  healing  rate  function,  co(q). 
For  large  doses  (ie.  long  times)  and  q  >  0,  the  exponential  term  is  small  and  can  be 
neglected,  making  ( I  h(q) !  2)t  independent  of  the  dose  Q.  From  the  experimental 
results  in  Figs.  4a  and  4b,  we  see  that  W  increases  more  like  Q  than  Q^2  and  that 
( I  h(q)  1 2>t  shows  a  significant  dose  dependence  at  large  q.  Thus,  although  linear 
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response  theory  qualitatively  accounts  for  a  number  of  the  observed  features,  it 
does  not  provide  us  with  a  quantitative  description  of  our  data. 

At  present,  no  non-linear  theory  exists  for  radiation  erosion,  but  a  general 
scaling  description^  has  been  developed  in  the  context  of  various  growth 
models.  According  to  scaling  theory: 

<  I  h(q)  1 2>t  oc  q-v  F(t  qz),  (4) 

with  F(x) xv/z  for  small  x,  and  F(x)  ~  constant  for  large  x.  This  would  predict  a 
power  law  decrease,  q*v  ,  in  <  I  h(q)  1 2)t  for  q  >  t*l^z.  In  Fig.  4a,  we  see  an 
approximate  power  law  dependence  in  the  correlation  function  at  large  q  for  the 
lower  ion  doses,  with  an  associated  exponent  of  order  -2.5  to  -2.9.  Under 
conditions  of  rotational  invariance,  the  exponent  z  is  related  to  v  by  z  =  2  -  a,  with 
2a  »  v  -  2,  resulting  in  z  -  1.6  - 1.8  and  v/z  -  1.5.  This  would  mean  that  for  small  q, 
<  I  h(q)  1 2)t  -  tv/z  should  increase  faster  than  linearly  with  increasing  time  (dose), 
while  the  correlation  length  at  crossover,  %  -  tl^z  ,  should  increase  more  slowly 
than  linearly.  As  we  have  seen,  both  predictions  are  valid  for  the  experiments  in 
this  study.  The  values  of  a  obtained  from  numerical  studies  in  three  dimensions 
are  model  dependent:  17  for  weak  non-linearity,  a  -  0.15  -  0.23  while  a  is  ~  4.0  for 
strong  coupling.  For  our  experiments,  we  find  a  to  be  -  0.2  -  0.4.  The  results 
embodied  in  Fig.  4  thus  appear  to  be  consistent  with  local  growth  models,  at  least 
at  lower  doses.  We  find  this  agreement  somewhat  surprising,  as  there  is  no  a 
priori  reason  to  expect  a  local  growth  model  to  correctly  describe  erosion  via 
sputtering.  16 

In  summary,  ion  bombarded  graphite  surfaces  evolve  a  rough  morphology 
characterized  by  the  divergence  of  the  correlation  length,  as  predicted  by  linear 
response  theory.  However,  the  experimental  height-height  correlation  function  is 
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not  quantitatively  consistent  with  a  linear  response  argument,  but  does  have 
features  in  common  with  the  scaling  theory  for  sputter  growth. 

This  work  was  supported  by  the  National  Science  Foundation  through 
Grant  DMR  8922027  and  by  the  Office  of  Naval  Research.  We  would  like  to  thank 
L.  Golubovic  and  R.P.U.  Karunasiri  for  helpful  discussions  and  E.  Snyder  for 
assistance  with  the  computer  software. 


8 


References 


1.  J.  L.  Vossen  and  W.  Kern,  Thin-Film  Processes  (Academic  Press,  New 
York:  1978) 

2.  H.  J.  Mathieu,  in  Thin  Film  and  Depth  Profile  Analysis,  H.  Oeschner,  ed. 
(Springer-Verlag,  Berlin:  1984),  pp.  39-58. 

3.  G.  Carter,  in  Erosion  and  Growth  of  Solids  Stimulated  by  Atom  and  Ion 
Beams,  G.  Kiriakidis,  G.  Carter,  and  J.  L.  Whitton,  eds.  (Martinus  Nijhoff, 
Hingham,  MA:  1986),  pp.  70-97. 

4.  I.  J.  Hodginson  and  P.  W.  Wilson,  CRC  Reviews  15,  27  (1988). 

5.  See,  for  example,  G.  Carter  and  M.  J.  Nobes,  Earth  Surface  Processes  5,  131 
(1980). 

6.  I.  V.  Katardjiev,  J.  Vac.  Sri.  Technol.  A  6  (4),  2434  (1988). 

7.  J.  A.  Kubby  and  B.  M.  Siegel,  in  Erosion  and  Growth  of  Solids  Stimulated  by 
Atom  and  Ion  Beams,  G.  Kiriakidis,  G.  Carter,  and  J.  L.  Whitton,  eds. 
(Martinus  Nijhoff,  Hingham,  MA:  1986),  pp.  444-452. 

8.  See,  for  example,  D.  M.  Mattox,  J.  Vac.  Sri.  Technol.  A  7,  1105  (1989). 

9.  G.  Carter,  B.  Navinsek,  and  J.  L.  Whitton,  in  Sputtering  by  Particle 
Bombardment,  Vol.  II,  R.  Behrisch,  ed.  (Springer-Verlag,  Berlin:  1983),  pp. 
231-266. 

10.  G.  K.  Wehner,  J.  Vac.  Sri.  Technol.  A  3,  1821  (1985). 

11.  M.  Kardar,  G.  Parisi  and  Y.  C.  Zhang,  Phys.  Rev.  Lett.  56,  889  (1986).  For  a 
review,  see  T.  Vicsek,  Fractal  Growth  Phenomena  (World  Scientific, 
Singapore,  1989). 

12.  F.  Family  and  T.  Vicsek,  J.  Phys.  A  18,  L75  (1985). 

13.  P.K.  Hansxna  and  J.  Tersoff,  J.  Appl.  Phys.  61  (2),  R1  (1987). 


9 


14.  For  earlier  measurements  of  <  I  h(q)  1 2)  by  microscopy,  see  G.  Rasigni  et.  al.. 
Surface  Science  162,  985  (1985). 

15.  J.  D.  Weeks,  in  Ordering  in  Strongly  Fluctuating  Condensed  Matter  Systems 
(Plenum,  New  York,  1980),  pp.  293. 

16.  G.  S.  Bales,  et  al. ,  Science  249,  264  (1990). 

17.  H.  Yan,  D.  Kessler  and  L.  M.  Sander,  Phys.  Rev.  Lett.  64,  926  (1990). 


10 


» 


Figure  Captious 

Fig.  1  Constant  current  STM  topograph  of  a  graphite  surface  after  sputtering 
with  a  flux  J  of  6.9  x  10*®  ions/cm^  sec  and  an  ion  dose  Q  of  10*® 
ions/cm^  at  room  temperature.  The  X  and  Y  dimensions  are  2400  A 
and  the  total  Z  dimension  is  18.6  A. 

Fig.  2  STM  topograph  with  the  same  parameters  as  that  in  Fig.  1,  except  that 
Q  =  1017  ions/cm^  and  the  total  Z  scale  is  27.9  A. 

Fig.  3  STM  topograph  with  the  same  parameters  as  that  in  Fig.  1,  except  that 
Q  =  10*®  ions/cm^  and  the  total  Z  scale  is  231.5  A. 

Fig.  4  Height-height  correlation  functions: 

(a)  (  I  h(q)  I  ^  >  for  the  data  in  Figs.  1-3,  with  Q  =  10*®  ions/cm^  ( O  ), 

10*7  ions/cm^  ( ■ ),  and  10*®  ions/cm^  (O).  The  crossover 
wavevector  q0  =  1/2;  is  indicated  for  each  curve  and  a  l/q2-7 
dependence  is  shown  for  comparison  in  the  large  q  regime.  The 
error  bars  indicate  one  standard  deviation  in  ( I  h(q)  I  ^ ) .  The 
uncertainty  decreases  as  q  increases  and  at  large  q  the  error  bars 
are  within  the  size  of  the  symbols  for  the  data  points. 

(b)  ( I  h(q)  I  ^ )  for  topographs  with  the  same  parameters  as  in  Fig.  4a, 
except  that  the  flux  is  higher:  J  =  3.5  x  10  ions/cm  sec.  Note  that 

n 

there  is  a  significant  flux  dependence,  with  ( I  h(q)  I  )  larger  (the 
surfaces  are  rougher)  for  the  higher  flux. 

(c)  Temperature  dependence  of  ( I  h(q)  I  ^ )  with  J  s  3.5  x  10*4  ions/cm^- 
sec,  Q  =  10*7  ions/cm^  and  T  =  300K  (□),  T  *  600K  (•)  and  T  *  900K  ( ♦ ). 
A  l/q4  dependence  is  shown  for  comparison  in  the  large  q  regime. 

Note  that  ( I  h(q)  I  ^ )  decreases  (the  surfaces  are  smoother)  with 
increasing  T. 
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